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The accelerated expansion of the Universe measured by high redshift Type la Supernovae ob- 
servations is explained using the non-equilibrium dynamics of naturally soft boson fields. Spinodal 
; ' instabilities inevitably present in such systems yield a very attractive mechanism for arriving at 

^ , the required equation of state at late times, while satisfying all the known constraints on models of 

■ quintessence. 

^ . One of the most startling developments in observational cosmology is the mounting evidence for the acceleration 
of the expansion rate of the Universe Coupled with cluster abundance and CMB observations, these data can be 
. interpreted as evidence for a cosmological constant A which contributes an amount 17a — 0.7 to the critical energy 
density while matter contributes fimattcr — 0.3, leading to a flat FRW cosmology (see [D and references therein). 

While introducing a cosmological constant may be a cosmologically sound explanation of the observations, it is a 
worrisome thing to do indeed from the particle physics point of view. It is hard enough to try to explain a vanishing 
cosmological constant, given the various contributions from quantum zero point energies, as well as from the classical 
theory H, but at least one could envisage either a symmetry argument (such as supersymmetry, if it were unbroken) 
or a dynamical approach (such as the ill-fated wormhole approach Q) that could do the job. It is much more difficult 
I ■ to see how cancellations between all possible contributions would give rise to a non-zero remnant of order 10~*^ GeV* 
y—{ ' which is extremely small compared to Mp^^^^^ or Afgugy, the "natural" values expected in a theory with gravity or 
, one with a supersymmetry breaking scale Msusy- Even from the cosmological perspective, a cosmological constant 
begs the question of why its effects are dominating now, as opposed to any time prior to today, especially given its 
[ different redshifting properties compared to matter or radiation energy density. 

• These fine tuning problems can at least be partially alleviated if instead of using a constant energy density to drive 
rS 1' the accelerated expansion, a dynamically varying one were used instead. This is the idea behind Quintessence models 
JL . A scalar field whose equation of state violated the strong energy condition (i.e. with p + 3P < ) during its 

Q ' evolution ||] would serve just as well as a cosmological constant in terms of explaining the data, as long as its equation 
(~| , of state satisfied the various known constraints for such theories [Q. However, an arbitrary scalar field whose energy 
density dominates the expansion rate is not sufficient to get out of all fine tuning problems; in particular, for a field 
of sufficiently small mass that it would only start evolving towards its minimum relatively recently (i.e. masses of 
order the inverse Hubble radius), the ratio of matter energy density pm to field energy density would be need to be 
?-j ' incredibly fine tuned at early times so as to have pm/ p<t> ~ 1 today. The quintessence approach uses so-called tracker 
fields Q that have potentials that drive the field to attractor configurations that have a a fixed value of Pm/Prj>- Thus, 
for these models, regardless of the initial conditions, the intermediate time value of Pm/p<p will always be the same. 
The only fine tuning required in these models is the timing of the deviation of from the tracking solutions to allow 
it to satisfy the condition p^ + SP^ < today. 

What we propose is a working alternative to the idea of tracker fields without the usual fine tuning problems. 
Recent work on the non-equilibrium dynamics of quantum fields P,p^ has shown that under certain circumstances 
the back reaction of quantum fluctuations can have a great influence on the evolution of the quantum expectation 
value of a field ||ll|,|l^, to the extent that using the classical equations of motion can grossly misrepresent the actual 
dynamics of the system. What we will show below is that we can make use of this modified dynamical behavior to 
construct models that might allow for a more natural setting for a late time cosmological constant. 

The class of models we consider are those using pseudo-Nambu-Goldstone bosons (PNGBs) to construct theories 
with naturally light scalars Such models have been used for late time phase transitions as well as to give 
rise to a cosmological "constant" that eventually relaxed to zero not unlike what we want to do here. However, 
our take on these models will be significantly different from that of ref. |p^ . 

We can write the required energy density as pDark Energy ~ (lO~^ eV) , which is suggestive of a light neutrino mass 



X 



scale [g^ 17|. There is a way to construct models of scalar fields coupled to neutrinos where the scalar field potential 
naturally (in the technical sense of t'Hooft [^) incorporates the small mass scale m'^. 
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Consider a Lagrangian containing a Yukawa coupling of the form ||T^ : 



i--Yuk = 2^ ymo + eexpj + j j VjlVjr + hj 



(1) 



The scale / is the scale at which a the global symmetry that gives rise to the Nambu-Goldstone mode $ is spontaneously 
broken. The Lagrangian /lyuk is to be thought of as part of the low-energy effective theory of $ coupled to neutrinos 
at energies below /. 

The term proportional to e could be obtained by a coupling to a Higgs field x that acquires an expectation value 
(x) = //V2expij. Note that in the absence of mg this Yukawa term possesses a continuous chiral U [V) symmetry. 
The term proportional to mg breaks this symmetry explicitly to a residual discrete Zj^ symmetry given by: 

27r/ 
N 

This interaction can generate an effective potential for the Nambu-Goldstone mode $ which must vanish in the limit 
that mo —f which is equivalent to the vanishing of the neutrino masses. Since $ is an angular degree of freedom, it 
should not be a surprise that the effective potential is periodic and of the form 

V ($) = [l + cos . (3) 

Here M should be associated with a light neutrino mass ruy ^ 10"'^ eV. 

Here, we have followed the working hypothesis of ||l5| which states that the effective vacuum energy will be dominated 
by the heaviest fields still evolving towards their true minimum. We assume that the super light PGNB field $, with 
associated mass of order ^ rn^/ f , will be the last field still rolling down its potential. Thus we have chosen by hand 
the constant in eq.(||) so that when <i> reaches the minimum it will have zero cosmological constant associated with 
it. This choice is essentially a choice of the zero of energy at asymptotically late times. 

The finite temperature behavior associated with these models is extremely interesting [|l9|. For TV > 3 the <f> 
dependent part of the potential can be written as 

c (T)m4 cos— , (4) 

where c (T) vanishes at high temperature T. Thus the high temperature phase of the theory has a non-linearly realized 
U {\) symmetry where the $ potential becomes exactly flat with value Af*. Since M ^ ra^ this cosmological constant 
contribution will have no effect during nucleosynthesis and through the matter dominated phase until T ~ M . At 
this time c(T) reaches its asymptotic value of unity and we have the potential in eq.(H). For iV = 2, c (T) changes 
sign continously, passing through zero at the critical temperature such that the high temperature minima become the 
low temperature maxima and vice- versa. There is a Z2 symmetry in both the low and high temperature phases. 

The potential in eq.(H) has regions of spinodal instability, i.e. where the effective mass squared is negative. These 
occur when cos A^$// > 0. If <i> is in this region, modes of sufficiently small comoving wavenumber follow an equation 
of motion that at least for early times is that of an inverted harmonic oscillator. This instability will then drive the 
non-perturbative growth of quantum fluctuations until they reach the spinodal line where cos iV$ / f ~ Q . Since 
the quantum fluctuations grow non-pertrubatively large, we have to resum perturbation theory to regain sensible 



behavior and this is done by the Hartree truncation |12|. The prescription is to first expand $ around its (time 
dependent) expectation value ($ {x, t)) = 4> (t) as 

$(f,t) = 0(t)+7A(f,i), (5) 

where the tadpole condition (x, t)) — gives the equation of motion for (j) (t). The Hartree approximation involves 
inserting eq.(|^) into eq.(^), expanding the cosines and sines, and then making the following replacements | [T^ : 

cos-^j^l jexp-^^, exp-^y^. (6) 

The equations for the field (j), and the fluctuation modes fk coupled to the scale factor a (t) are 

0-h3-0 ^exp ^75-^ sm -r = 0' v) 

a f 2/^ / 



fk + S-Jk + 75— exp - 
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with 



(8) 



The effective Friedmann equation for the scale factor is obtained by use of semiclassical gravity, i.e. by using (T^i/) 
to source the Einstein equations: 



3M2 



Pra{t) + 2 + 2 + ^ ( 7 + ^^M 1 + cos(Af<^//) exp AL^ 



with 



Pm{t) = Prn{ti) 



aHt) 



(9) 



(10) 



being the matter density and ti being the time at which the PGNB field begins its evolution. 

The interesting feature of the above equations of motion is the appearance of terms involving exp - 
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These 

multiply terms in the potential and its various derivatives that contain the non-trivial <j) dependence. What we expect 
to have happen is that as the spinodally unstable modes grow, they will force (V'^) to grow as well. This in turn will 
rapidly drive the exponential terms to zero, leaving a term proportional to M'^ in the Friedmann equations, which 
will act as a cosmological constant at late times. 

If we consider the iV > 3 models, then at temperatures larger than Tdit ~ M the potential is just given by 
and is swamped by both the matter and radiation contributions to the energy density. Since the potential is flat, 
we expect that the zero mode is equally likely to attain any value between and 27r and in particular, we expect 
a probability of order ^1/2 for the initial value to lie above the spinodal line. If there was an inflationary period 
before this phase transition, we expect that the zero mode will take on the same value throughout the region that 
will become the observable universe today. 

As the temperature decreases, the non-trivial parts of the potential turn on and the zero mode begins its evolution 
towards the minimum once the Universe is old enough, i.e. H (Troii) ~ ~ / f . At the same time, if the zero 
mode started above the spinodal line, the fluctuations begin their spinodal growth. Whether the spinodal instabilities 
have any cosmological effect depends crucially on a comparison of time scales, the first being between the time 
it takes the zero mode to reach the spinodal point at 0spinodai// = 7r/2iV under the purely classical evolution (i.e. 
neglecting the fiuctuations), and the time ^spinodal it takes for the fluctuations to sample the minima of the tree-level 
potential, so that N'^{'4>^) I P ~ C (!)■ Since the growth of instabilities will stop at times later than t,, if spinodal 
instabilities are to be at all relevant to the evolution of 0, we need ^spinodal ^ t*- By looking at the equations of 
motion we can argue that 



□ dal 
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In 



2M2 iV2(?/'2) (ii) 2H, 



(11) 



where ti is the time at which the zero mode starts to roll and Hi is the Hubble parameter at this time. Furthermore 
the early time behavior of the equations of motion gives us 



M2 ^ N(j){U) 



+ 



2H,, 



(12) 



Comparing eqs.( [ri| , [r2[) , we see that to have the spinodal instabilities be significant we need (jP'iti) <^ (ip^) {ti). 

The other condition that needs to be met is that there should be sufficient time for the spinodal instabilities to 
dominate the zero mode evolution before today. This will ensure that the expansion of the Universe will be driven 
by the remnant cosmological constant M'^ at the times relevant to the SNIa observations. For large enough initial 
fluctuations we can make the spinodal time as early as we need. 

What sets the scale of the initial fluctuations? If we assume a previous inflationary phase, we can treat the PGNB 
as a minimally coupled massless field and the standard inflationary results should apply EG]. The initial conditions 
for the mode functions are then given by: 



fk iu) 



-Hds for K < k < Hi. 



(13) 
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where k is an infrared cutoff corresponding to horizon size during the De Sitter phase and Hds is the Hubble parameter 
during inflation. The short wavelength modes (fc > Hi) have their conformal vacuum initial conditions. With these 
initial conditions 
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(A^c-toids - 60) . (14) 



and for Hds - 10^^ GeV, A^c-foids - 60 ~ 1G^ and / ~ lO^^ GeV, we would only need that (j}{U)/f < 0.5. These are 
not outlandish parameter values and we see that very little fine tuning is required. In fact, there is a great deal of 
freedom in choosing the values of these parameters, the only requirement being that the ratio of parameters appearing 



in (14) not be so small that the required value of (j){ti)/ f is overly restricted. 

In the figures below we use these parameters as well as M ~ 5.5 x 10""^ eV corresponding to neutrino masses, 
beginning the evolution at a redshift 1 + z = 1200. In Fig. |l| we plot the numerical evolution of the zero mode and 
of the growth of the fluctuations, while Fig. || shows the equation of state of the PNGB field and the total equation 
of state including PNGB and matter components as a function of redshift. What we quickly infer from these graphs 
is that the evolution of the Universe becomes dominated by the remnant cosmological term, leading to an evolution 
toward a late time equation of state w = P/ p ~ — 1. The equation of state today is seen to be w ~ —0.7 and indicates 
a matter component fimattcr = 0.3 and a cosmological constant-like component f2pngb = 0.7. Because the PNGB 
component has an equation of state w = —1 by a redshift as high as z = 50, these results reproduce the best fit 
spatially flat cosmology of the SNIa data . 

One feature of this model is that the parameter M is directly related to the measured value of today's Hubble 
constant. We find 



which is to be compared to the observed 90% confidence range of Am^ from the Super Kamiokande contained events 
analysis ||l^ of 5 x lO^^eV'^ < Ato^ < 6 x 10^'^eV'^, and the more recent results of the up-down asymmetry analysis 
which indicates a range 1.5 x lO^'^eV^ < Am^ < 1.5 x 10~^eV^, whereas the small and large mixing angle MSW 



solutions of the solar neutrino problem |21 2^ yield a range of 5 x 10~^eV^ < Am^ < 4 x lO^'^eV^. 

There is no shortage of models to explain the accelerating expansion of the universe. However, most options are 
lacking in motivation and require significant fine tuning of initial conditions or the introduction of a fine tuned small 
scale into the fundamental Lagrangian. We too have a fine tuned scale: the neutrino mass. However, we can take 
solace in the fact that this fine tuning is related to a particle that can be found in the Particle Data Book |Q , with 
known mechanisms to produce the required value, and experiments dedicated to its measurement. 

The model itself is also relatively benign, not requiring invocations of String or M-theory to justify its potential. 
Chiral symmetry breaking leading to PNGB's is not unheard of in nature (pions do exist after all!), and should probably 
be expected in GUT or SUSY symmetry breaking phase transitions involving coupled scalars. This, together with 
the dynamical effects of backreaction allow the present model to be successful in explaining the data with only minor 
tuning of initial conditions. 
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FIG. 1. The evolution of the zero mode (f> and the root mean square flucutation {ijj'^)^^'^ as a function of the redshift 1 + z 
from recombination at 1 + z ~ 1200 to today at 1 + ^=1 and beyond. The parameters used for the simulation correspond to 
M = 5.5 X 10~^ eV, f/N = 10^^ GeV and (A^e-foids - 60)iyDs = 10^^ (GeV)^ For convenience we have used (A^c-foWs - 60) = 10, 
which is necessary to avoid numerical implementation of an infrared cutoff of order exp( — 10^). This has the effect of slightly 
underestimating {-ip^) as the evolution proceeds compared to the case of (A^e-foids — 60) — 10^, but will not change any qualitative 
features of the evolution. The initial conditions are N(j){ti)/f = 0.08 and (j}{ti) = 0. 
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FIG. 2. The total equation of state lutotai ~ Ppngb / {pmattcT + ^ngb) and the PNGB equation of state xOpngb ~ f'pngb/ppngb 
as a function of redshift. The parameters were chosen as in Fig. 
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